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We present a measurement of the inclusive jet cross section using the Run II cone algorithm and 
data collected by the DO experiment in pp coUisions at a center-of-mass energy ^/s =1.96 TeV, 
corresponding to an integrated luminosity of 0.70 fb~^. The jet energy cahbration and the method 
used to extract the inclusive jet cross section are described. We discuss the main uncertainties, which 
are dominated by the jet energy scale uncertainty. The results cover jet transverse momenta from 
50 GeV to 600 GeV with jet rapidities in the range —2.4 to 2.4 and are compared to predictions using 
recent proton parton distribution functions. Studies of correlations between systematic uncertainties 
in transverse momentum and rapidity are presented. 

PACS numbers: 13.87.Ce, 12.38.Qk 



I. INTRODUCTION AND MOTIVATION 



The measurement of the cross section for inclusive pro- 
duction of hadronic jets in hadron coUisions provides 
stringent tests of quantum chromodynamics (QCD). The 
inclusive jet cross section in pp collisions for jets with 
large momentum transverse to the beam axis (pt) is di- 
rectly sensitive to the strong coupling constant (as) [l| 
and the parton distribution functions (PDFs) of the pro- 
ton 0. At the Tevatron pp collider, data are divided 
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into two sets corresponding to Run I (1992 — 1996) and 
Run II (2002 — 2011). The increased pp center-of-mass 
energy between Run I {y/s — 1.8 TeV) and Run II (-y/s = 
1.96 TeV) leads to a significant increase in the cross sec- 
tion at large px - a factor of three at px ~ 550 GeV, 
as shown in Fig. [T] obtained using the next-to-leading 
order (NLO) QCD calculation as implemented in nlo- 
JETH — h 01 ■ This increases the sensitivity to potential 
new observations such as quark compositeness and ex- 
tra dimensions [4|. The integrated luminosity of the in- 
clusive jet cross section measurement discussed in this 
Article exceeds the Run I luminosity by more than a fac- 
tor of five, allowing for more stringent constraints on the 
PDFs. In Fig.[5]we show the different subprocesses that 
contribute to the inclusive jet cross section. In particular 
the gluon density can be further constrained using these 
data, since the gg and qg initial states contribute signifi- 
cantly to the cross section across almost the full px range 
of the measurement. The gluon distribution is still poorly 
known, especially for gluons carrying a large momentum 
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FIG. 1: (color online) Inclusive jet production cross section 
for central jets {\y\ < 0.4) for Run I and Run II energies 
at the Tevatron obtained using NLO QCD as implemented 
in nlojetH — h. The ratio of the two curves is shown in the 
bottom panel. We note an increase of the Run II cross section 
with respect to Run I of up to a factor 10 at highest jet pr- 

fraction x. In contrast, the quark PDFs are already well 
constrained by fixed target and electron-proton collider 
experiments [2]. 

In this Article, we report measurements by the DO col- 
laboration of the inclusive jet cross section in pp collisions 
at a center-of-mass energy of -^/s = 1.96 TeV. We give 
details of the analysis leading to the results published 
in Ref. 0] , with particular attention to the jet energy 
scale determination. The precision achieved for the jet 
energy scale in the DO experiment is unprecedented for 
any hadron collider experiment to date, and the meth- 
ods applied to reach this precision will be useful for future 
hadron collider experiments. 

The data sample, collected with the DO detector dur- 
ing 2004-2005 in Run II of the Fermilab Tevatron, corre- 
sponds to an integrated luminosity of£ = 0.70fb~^ 0. 
The cross section is presented in six bins of jet rapid- 
ity (y), extending to |y| = 2.4, as a function of jet px 
starting at pr = 50 GeV. The rapidity is related to the 
polar scattering angle 9 with respect to the beam axis by 
y = O.51n[(l+^cos6i)/(l-/3cos6')] with (3 = \p\/E. The 
measurement also extends the kinematic reach of earlier 
measurements of the inclusive jet cross section by the 
CDF and DO CoUaborations ^U^. 

This Article is organized as follows. After a brief de- 
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FIG. 2: (color online) Fractional contributions of the qq, qg 
and gg sub-processes to the inclusive jet cross section for cen- 
tral jets as a function of jet pr and of the fraction of the beam 
energy carried by the jet, xt- 

scription of the DO detector in Sec. [Ill we discuss the jet 
algorithm used in Run II in Sec. lIIIl Section lTVl describes 
the theoretical predictions for the inclusivejet cross sec- 
tion before the DO measurement in Ref. [5| . Section |V] 
gives an extensive description of the methods used to 
measure the jet energy scale and to determine the corre- 
sponding uncertainty. This is the leading uncertainty for 
the measurement of the inclusive jet cross section. Sec- 
tions |Vl]-|Xl] describe the jet triggers, event and jet selec- 
tion criteria, determination of the jet px resolution and 
the unfolding method. In Sections IXIII IXIVl we describe 
our results and compare them with predictions using re- 
cent PDF parameterizations. 



II. DETECTOR 

In this section, we briefly describe the Run II DO de- 
tector [l^l and the main components used in the mea- 
surement of the inclusive jet cross section. 



A. Calorimeter 

The calorimeter and the tracking detectors, used to 
measure the position of the interaction point, are the 
most important detector components used to measure 
the jet Pt- An accurate and stable energy response is re- 
quired for reliable measurements of the cross section for 
jet production. The calorimeter consists of the follow- 
ing subdetectors: the uranium/liquid argon calorimeter 
divided into a central (CC) and two end (EC) sections, 
the plastic scintillator inter-cryostat detector (ICD), and 
the massless gap (MG) detectors. Both the CC and ECs 
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FIG. 3: Side view of a quadrant of the DO calorimeters (CC, EC and ICR) showing the transverse and longitudinal segmen- 
tation [13]. The shading pattern indicates the cells for signal readout. The lines indicate the pseudorapidity intervals defined 
from the center of the detector. The CC covers the region \ri\ < 1.2 and the EC extends the coverage up to \ri\ ~4.2. The 
inter-cryostat detector is visible as a thin dark shaded tile between the cryostats, within 1.1 < \ri\ < 1.4, and the massless gap 
detectors are inside the cryostats, within 0.8 < \r]\ < 1.2 (in the CC) and 1.0 < \ri\ < 1.3 (in the EC). 



are segmented longitudinally into electromagnetic (EM) , 
fine hadronic (HAD), and coarse hadronic (CH) sections. 
A schematic view of the calorimeter showing its pro- 
jective tower geometry as a function of pseudorapidity 
rj = — lntan(6'/2), where 9 is the polar angle from the 
beamline, is given in Fig. [3l The choice of binning in 
the inclusive jet cross section measurement closely fol- 
lows the structure of the calorimeter: \ri\ < 0.8 is well- 
contained within the CC, 1.6 < \ri\ < 2.4 within the EC, 
whereas the more challenging inter-cryostat region (ICR) 
0.8 < 1 77 1 < 1.6 has energy sharing between the four sub- 
detectors. 



1. Central and end calorimeters 

The CC covers detector pseudorapidity I77I < 1.2, 
and the two ECs extend the range up to \r]\ = 4.2. 
Both the electromagnetic and fine hadronic calorime- 
ters are sampling calorimeters with an active medium 
of liquid argon and absorber plates of nearly-pure de- 
pleted uranium. Incoming particles traversing the ura- 
nium absorber plates initiate showers of secondary par- 
ticles that ionize the argon in the gaps between the ab- 
sorber plates. A high-voltage electric field collects the 
free electrons on resistively-coated copper pads that act 
as signal boards [l^, [ill • The outer part of the calorime- 
ter, the coarse hadronic section, uses copper in the CC 
and stainless steel in the EC for the absorber plates. The 
calorimeter is transversely segmented into cells in pseudo- 



rapidity and azimuthal angle of 0.1 x 0.1 (and 0.05 x 0.05 
in the third layer of the EM calorimeter) for |?7| < 3.2 to 
allow for more precise location of EM shower centroids. 
At I?/! > 3.2, the cell size grows to 0.2 or more for both 
Tj and the azimuthal angle 4>. These high pseudorapidi- 
ties are not used for the jet cross section measurements 
since the jet triggers are limited to \ri\ < 3.2. The total 
depth of the EM calorimeter is about 20 electromagnetic 
radiation lengths, and the combined thickness of the elec- 
tromagnetic and hadronic calorimeters is about 7 nuclear 
interaction lengths [l3{ . 



A typical calorimeter cell consists of an absorber plate 
and a liquid argon gap. The metal plate is grounded, 
while the resistive plate of the signal board located in 
the liquid argon gap is kept at a high voltage of ~2.0 kV. 
The drift time of the electrons across the typical 2.3 mm 
gap is 450 ns, longer than the separation between two 
subsequent Tevatron bunch crossings of 396 ns. To mini- 
mize the effect of pile-up from interactions from different 
bunch crossings, only two-thirds of the charge collected is 
used in the shaper circuits and then provided to baseline 
subtraction boards. To remove the baseline, the signal 
corresponding to a sampling occurring 396 ns earlier (the 
time between two bunch crossings) is subtracted. Only 
cells with a signal at least 2.5 times the standard devia- 
tion of the electronic noise after baseline subtraction are 
kept in nominal conditions of data taking. This defines 
the on-line zero-suppression mode of the calorimeter. 
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2. Inter-cryostat detector and massless gaps 

The regions between the CC and the ECs are in- 
strumented with the inter cryostat detector and mass- 
less gaps. The ICD and MG detectors provide energy 
measurement for the otherwise poorly instrumented re- 
gions located at roughly 0.8 < \ri\ < 1.4, where the 
depth of the passive material coming from cryostat walls, 
stiffening rings and cables varies rapidly with rapidity. 
The ICD relies on photomultipliers to record the sig- 
nals from plates of scintillating plastic and covers the 
region 1.1 < \ri\ < 1.4. The signal from the ICD is 
stretched in time to match that of the EM calorime- 
ter and augments the EM calorimetry that is absent in 
the region 1.2 < \?]\ < 1.35. The ICD is supplemented 
by the MG detectors that are placed inside the cryostat 
walls in the CC and the ECs from 0.8 < |?7| < 1.2 and 
1.0 < |?7| < 1.3, respectively. Unlike typical calorimeter 
cells, the massless gap detectors do not have absorber 
plates, but they sample the showers that develop in the 
cryostat walls, calorimeter support structures, and other 
calorimeter cells. 

In addition to the CC, ECs, and ICD, preshower detec- 
tors are located in the central and forward regions, but 
they are not used in this analysis. 



B. Tracking detectors 

The tracking detectors are not used directly in jet re- 
construction since the jet finding algorithms in DO use 
only energy deposits in the calorimeter towers. However, 
the tracking detectors are used to reconstruct the posi- 
tion of the primary vertex of the pp interaction, which is 
necessary to precisely measure the jet rapidity and trans- 
verse momentum. The position of the primary vertex is 
typically distributed as a 20 cm-wide Gaussian distribu- 
tion along the beamline direction around the nominal 
interaction point of (x,y,z) = (0,0,0) located in the 
center of the detector. In the detector description and 
data analysis, we use a right-handed coordinate system 
in which the z-axis is along the proton direction and the 
y-axis is upward. The inner tracking system, consisting 
of the silicon microstrip tracker, provides a 35 /um vertex 
resolution along the beam line and 15 /im resolution in 
the r — (f> plane for tracks with a minimum px of 10 GeV 
at rj = 0. The outer tracking system, consisting of the 
central fiber tracker, uses scintillating fiber technology 
to complement the silicon tracker. Both detectors are lo- 
cated in the 2 T magnetic field of the superconducting 
solenoidal magnet to allow measurements of the momen- 
tum of charged particles. 



C. Muon detector 

The muon detector is composed of a combination of 
proportional drift tubes in the central region {\ri\ < 1.0), 



and smaller, faster mini drift tubes in the forward region 
(1.0 < \ri\ < 2.0). Both are separated in three layers (A, 
B, C). Toroidal magnets are located between the A and 
B layers of the muon detector in the central and forward 
regions to allow reconstruction of the muon momentum. 
The muon system is not used directly in our analysis (we 
do not correct for muons in jets), but very high energy 
jets can leak outside the calorimeter and show some hits 
in the A layer. We do not include these hits in jet re- 
construction, but instead correct the jet cross sections 
for asymmetries introduced in the jet energy resolution 
(described in Scc UXCl) . 



D. Luminosity detector 

The luminosity monitor (LM) is constructed of scin- 
tillating tiles on both sides of the interaction point that 
detect the particles coming from inelastic collisions. The 
luminosity C is determined from the average number of 
observed interactions -/Vlm using the formula 



LM 



C^LM 



(1) 



where / is the pp bunch crossing frequency, and ctlm 
is the effective cross section for inelastic collisions mea- 
sured by the LM that takes into account event losses due 
to inefficiencies and geometric acceptance [tJ. In prac- 
tice, A'lm is calculated by inverting the expression for 
the Poisson probability of observing zero LM hits in ei- 
ther of the two arrays 



P(0) 



-CTLM^// y, ^2e~°'ss£/(2/) _ g-<Tss£//^ ^ (2) 



The right-most term of Eq. [5] accounts for the possibility 
of producing double-sided LM hits from a combination of 
single-sided (SS) LM hits, where ass is the effective cross 
section for only one of the arrays to show hits. The un- 
certainty on the luminosity determination is estimated 
to be 6.1% Q. This uncertainty is dominated by the 
5.4% uncertainty coming from the determination of ctlm, 
roughly half of which is due to acceptance and efficiency 
of the LM detectors with the remainder due to the un- 
certainty in the total inelastic cross section at 1.96 TeV 
described in [7l.[Tl|. 



III. JET RECONSTRUCTION 

Jets are reconstructed using the Run II midpoint cone 
algorithm p^ . which is an iterative cone algorithm that 
considers energy deposits as four- vectors to construct the 
jet four-momentum. The same algorithm is used with dif- 
ferent inputs in data and Monte Carlo (MC). It is used 
to build jets from energy deposits in the calorimeter in 
data or in fully simulated MC events, out of stable par- 
ticles in simulation, and out of partons produced either 
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in a parton shower simulation or from a next-to-leading 
order theoretical calculation. 

In data and in MC events processed through a simula- 
tion of the response of the DO detector, the first step is to 
define the seeds for jet reconstruction. Pseudoprojective 
towers, as illustrated in Fig. [U are built by adding the 
4-momenta of the calorimeter cells. The 4-momentum 
associated with the energy deposit in each cell of the 
calorimeter is computed using the direction defined by 
the reconstructed pp interaction vertex and the center of 
the cell and assuming E — \p\. All non-zero-suppressed 
cells are used in jet reconstruction. Calorimeter towers 
are ordered in decreasing transverse momentum and are 
used as seeds to form preclusters using a simple cone al- 
gorithm of radius 0.3 in (77, (jj) plane, starting with the 
tower having the highest px and then descending the list 
until no towers remain above a minimum threshold of 
PT > 500 MeV. All towers added to a precluster are re- 
moved from the list, avoiding overlaps between preclus- 
ters. Preclusters with pr > 1 GeV are used as seeds for 
the jet clustering algorithm. The goal of preclustering in 
data is to reduce the number of seeds and the computing 
time to reconstruct jets. As verified by MC studies [l5| . 
the low value of the pt threshold on the jet seeds ensures 
that there are no significant variations in the jet observ- 
ables for the pr range considered in this measurement 
(pT > 50 GeV). 

The seeds - preclusters in data and in MC events 
processed through a simulation of the response of the 
DO detector, or stable particles in MC, or partons from 
NLO calculation - are used as center points for proto- 
jets. All calorimeter towers, particles or partons within 
AR = ^(Ay)2 + (A(/.)2 < where i?cono = 0.7, 

are added to the proto-jet. The four-momentum of the 
proto-jet is the sum of the four-momenta of all included 
calorimeter towers, particles or partons. The direction of 
the resulting four- vector is used as the center point for a 
new cone. When the proto-jet four-momentum does not 
coincide with the cone axis, the procedure is repeated 
using the new axis as the center point until a stable 
solution is found. The maximum number of iterations 
is 50 and the solution is considered to be stable if the 
difference in AR between two iterations is smaller than 
0.001. In the rare cases of bistable solutions the last iter- 
ation is retained. Any protojets falling below a threshold, 
^'rjet ^ ^'Tmiu' '^itli Prmin ~ ^ GeV, are discarded. 

The presence of a threshold requirement on the cluster 
seeds introduces a dependency on infrared and collinear 
radiation. In order to reduce the sensitivity to soft radia- 
tion, pr-weighted mid-points between pairs of proto-jets 
are used as additional seeds if the distance between pairs, 
AR in the (y, (j>) plane to the proto-jet, is between 0.7 
and 1.4. The list of stable proto-jets obtained from this 
procedure may contain many overlapping and identical 
jet candidates. To resolve these ambiguities the proto- 
jets are sorted in order of decreasing px and processed 
through a split-and-merge procedure to remove overlaps. 
If two proto-jets have overlapping cones, they are merged 



if the overlap region contains more than 50% of the trans- 
verse momentum of the lower px jet. Otherwise, the jets 
are split with calorimeter cells or particles in the over- 
lap region being assigned to the nearest jet in (y, 0). In 
both cases, the jet four-momenta are recomputed after 
this reassignment. In case of multiple overlaps, the al- 
gorithm always starts with the highest px proto-jet to 
redistribute the shared towers. As mentioned above, the 
jet four-momentum is computed as the sum of the four- 
momenta of the (massless) calorimeter energy deposits 
included in the jet, and consequently the calorimeter jets 
are massive by construction if the jet cone contains cells 
with different locations in the {rj, 0) plane. The variables 
used to characterize the jets are the jet pr and y. The 
split-and-merge procedure may modify the cone axis and 
jet four-momentum for the final jets, and include towers 
outside the initial 0.7 cone. 



IV. THEORETICAL PREDICTIONS 

In this section, we describe how we compute the pre- 
dictions of the inclusive jet cross sections that are later 
compared to our measurements. 

1. Jet cross section at NLO 

We use the program FastNLO [l^, which is based 
on the matrix elements implemented in NLOJET-I — h, to 
calculate the inclusive cross sections to next-to-leading 
order precision and to evaluate the effects of the choice of 
proton PDFs, such as CTEQ6 or MRST2004 [Ulii, in 
a computationally efficient manner. Perturbative QCD 
(PQCD) requires the specification of the renormalization 
scale /ii? and the factorization scale fip. Typical choices 
set both pfi — i^ip to the pT of each of the individual jets, 
with half and twice this scale used to estimate the the- 
oretical scale uncertainty. The uncertainty on the NLO 
prediction of the inclusive jet cross section due to the 
choice of renormalization and factorization scales is given 
in Fig. Hand is about 10-20%. 

2. Parton distribution functions 

A discussion of the different PDFs and methods to re- 
duce their uncertainties using various measurements at 
the Tevatron and the LHC can be found in reference 
In this paper, we bricfiy describe the PDFs used in the 
comparison between the measurements and the theoret- 
ical predictions. 

One of the PDF sets used in this analysis is pro- 
vided by the CTEQ Collaboration. This most re- 
cent global fit from the CTEQ Collaboration performed 
prior to the measurement described in this paper, called 
CTEQ6.5M [3, utilizes DO and CDF Run I measure- 
ments, as well as the most recent deep inelastic scatter- 
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FIG. 4: (color online) Uncertainty on the inclusive jet cross section due to the choice of the renormalization and factorization 
scales fiR and /if in the NLO QCD calculation using NLOJET+ + . 



ing (DIS) data from the HERA collider at DESY and 
existing fixed target DIS and Drell-Yan data. The cen- 
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tral prediction of the CTEQ6.5M PDF is supplemented 
with the provision of 20 eigenvector basis PDF sets to 
estimate the PDF uncertainty, representing independent 
variations of the PDFs virithin the 90% C.L. of the data 
sets used in the fit. 

Another widely used PDF parameterization is pro- 
vided by the MRST Collaboration 1^. Our measure- 
ments are compared to the MRST2004 parameterization, 
which does not include our results. A third PDF param- 
eterization is MSTW2008 JJj which uses our results. 
The differences with respect to CTEQ6.5M are mainly 
in the description of the gluons at high-x and are within 
the CTEQ6.5M uncertainty band, as shown in Fig. [S] 
We also note that the uncertainty on the gluon density 
calculated by the CTEQ6.5M parameterization is larger 
than 40% for x > 0.5 and squared four-momentum trans- 
fer — 500^ GeV. Comparisons between our data and 
NLO calculations using these and other PDF parameter- 
izations are given in Sec. IXIIl 



V. JET ENERGY SCALE MEASUREMENT 

In this section we describe the method used to obtain 
the jet energy scale (JES) applied in the measurement of 
the inclusive jet cross section as a function of jet px- To 
compare the theoretical predictions to data, both need to 
be corrected to a common reference-level, chosen here to 
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be the "particle-level jets." We correct the calorimeter jet 
energies to the particle level, and apply non-perturbative 
corrections (hadronization and underlying event) to theo- 
retical NLO cross sections to move from the parton to the 
particle level. Particle jets (201 are clustered from stable 
particles after fragmentation, including particles from the 
true underlying event, but excluding undetected energy 
from muons and neutrinos. The JES procedure provides 
a correction factor that translates on average the energy 
of jets measured in the calorimeter to the energy of the 
corresponding particle jets. The jet energy scale is de- 
termined from data acquired during the same running 
period as used in the measurement of the inclusive jet 
cross section. 

The main effects that need to be considered when cor- 
recting jet energies from the calorimeter measurement 
Emcas to the particle level -Eparticic are the offset energy 
(O), calorimeter response (i?), and detector showering 
(S). These corrections can be expressed as a simple for- 
mula 



The offset energy O originates from electronics noise, 
calorimeter noise from uranium decays, residual energy 
from previous bunch crossings ("pile- up"), and energy 
from multiple pp collisions during a bunch crossing. The 
underlying event energy corresponding to multiple par- 
ton interactions in a single pp collision is not considered 
as part of the offset energy since it is included in the jet 
energy at the particle level. This also avoids correcting 
the data with model dependent offset corrections. The 
calorimeter response R is the average fraction of the en- 
ergy measured in the calorimeter for the particles inside 
the jet cone. The detector showering is the net flow of 
energy in and out of the jet cone due to detector effects, 
such as the magnetic field, scattering from passive mate- 
rial, and shower development in the calorimeter. The cor- 
rection S is defined as the ratio of the response-corrected 
calorimeter jet energy, in the absence of offset, and the 
particle jet energy. The correction does not include the 
effects of real QCD emissions, which arise from partons 
that shower outside the jet cone. We discuss each correc- 
tion in turn below. 



The typical value of the NP offset in a cone, R = 0.7, is 
0.2 GeV in the CC and ECs and 0.5 GeV in the ICR for 
the instantaneous luminosities considered in this analy- 
sis. 

The MI offset is the energy deposited by additional 
collisions during the bunch crossing. The value of the MI 
offset increases linearly with the number of additional 
interactions, which is characterized by the number of re- 
constructed pp interaction vertices in a given event. A 
typical value of MI is of the order of 0.5 GeV in the CC 
per additional interaction. 

The offset energies are measured directly from data us- 
ing "zero bias" and "minimum bias" data collected at a 
constant rate of about 0.5 Hz during data taking. The 
only requirement for zero bias events is coincident timing 
with the beam crossing; minimum bias events addition- 
ally require energy depositions above thresholds in coin- 
cidence in the two luminosity monitors, indicating that 
an inelastic collision took place. The offset is estimated 
from the average energy density in all calorimeter towers 
within detector rings of fixed pseudorapidity. The offset 
energy for a given jet cone is then calculated by sum- 
ming the average offset in towers within the cone radius 
around the jet center. The NP offset energy is measured 
using zero bias data with a veto on the luminosity mon- 
itor (no interaction occurred), and the MI energy for a 
given number N of interactions is the difference in the 
energy in minimum bias events with (N + 1) vertices and 
with a single vertex. 

The offset energy for different numbers of pp interac- 
tions (measured by the number of reconstructed vertices) 
is displayed in Fig. [6] and is found to depend linearly on 
the number of interactions within a 5% uncertainty. The 
average vertex multiplicity in the sample used to measure 
the inclusive jet pr cross section is ~ 1.5 - 2.0, hence the 
average offset correction to jet pr is ~0.5 GeV in the CC 
and EC and ~0.7 GeV in the ICR. The uncertainties on 
the offset corrections are of the order of 1% of the overall 
energy correction at low jet pt and are negligible for jet 
PT above ~100 GeV. They are significantly smaller than 
the total jet energy scale uncertainties. 

B. Determination of the jet energy response 



A. Determination of the offset energy 

The offset energy consists of the energy in the jet that 
is not related to the primary pp collision (hard scatter and 
underlying event). The offset energy is divided into two 
distinct categories, noise and pile- up (NP), and multiple 
pp interactions (MI). The noise component corresponds 
to the contributions of calorimeter and electronics noise, 
as well as the decay of the uranium nuclei in the calorime- 
ter. The pile-up energy corresponds to the energy left in 
the calorimeter from previous or next collisions because 
of the long integration time of the calorimeter electronics. 



The jet energy response, i?, can be factorized into two 
parts R — Rcc{E) ■ F^{r],E). The Rcc term uses the pT 
balance between the 7 and the jet in 7-l-jet events with 
a high (photon) purity in the CC region to determine an 
absolute response correction, while the second term i*"^ 
normalizes the response of the calorimeter as a function 
of jet pseudorapidity. 

1. Jet response in the CC 

The missing transverse energy {^t) projection fraction 
(MPF) method [2l[ is applied in 7 -I- jet events to mea- 
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FIG. 6: (color online) Offset corrections as a function of the 
jet pseudorapidity in the detector (without taking into ac- 
count the vertex position) for different numbers of recon- 
structed primary vertices A'pv- The special case A^pv = 1 
includes only the noise contribution to the offset. 



sure the response for jet energies in the CC region. Use 
of the MPF reduces the sensitivity of the measurement 
to showering and additional unreconstructed jets. We 
project the vector sum of all calorimeter tower energies 
transverse to the beam (including those of the photon), 
which equals the opposite of the l^x in the event, onto the 
photon transverse momentum vector pt,^- At the parti- 
cle level, the photon is balanced against the hadronic 
recoil, pT,-y + PT,had — 0, where pT,-y and pr.had are the 
transverse momentum of the photon and the hadronic 
recoil system, respectively. The measured jet pr will be 
affected by the energy response of the calorimeter causing 
an imbalance in the jet and photon transverse momenta, 
resulting in a non-zero l^xi 



Rem ■ PT,-f + Rha.d ' PT,had = — -^Ti 



(4) 



where Rem and i?had are the electromagnetic and 
hadronic calorimeter responses, respectively. 

The MPF method necessitates a precise energy calibra- 
tion for electrons and photons. The electron energy scale 
is determined from data using Z — >■ e+e" decays [2^ . 
MC simulations tuned to reproduce the response for elec- 
trons in data are used to derive the response difference 
between photons and electrons. The leading uncertainty 
in this simulation is caused by limited knowledge of the 
number of radiation lengths of material in front of the 
calorimeter. 

Using the corrected photon energy scale (i?om — 1), 
-Rhad is determined after projecting all terms in Eq. Ulon 
the photon px unit vector fij. In the MPF method, the 



ad 



\PT,' 



(5) 



where we use Ipr,- 



Pr.had- When the jet is re- 



quired to be back-to-back with the photon (difference in 
azimuthal angle larger than 2.9 radians) and no addi- 
tional jets are allowed in events with a single pp interac- 
tion, the hadronic recoil response i?had can be identified 
with the jet response i?jct- The impact of the proton 
remnants is small on average. The jet energy response 
depends on the particle jet energy and the results are usu- 
ally binned in jet pt- However, the measured jet energy 
has poor resolution and can lead to a large bias in the 
measurement of the response. To avoid this resolution 
bias, the jet energy response is measured as a function of 
the estimator 



E' = pT.-y ■ coshr^jc 



(6) 



E' is strongly correlated to the particle level jet energy 
and has a better resolution than the measured jet energy. 
We parameterize all corrections as a function of E' and 
map back to the measured jet energy i^moas on a jet-by- 
jet basis by inverting the equation 



-E^mcas — O — i?had(-B')5phvs(-B')i?', 



(7) 



where O is the offset contribution, i?had(£'') contains 
all jet energy corrections back to particle level, and 
Sphys{E') = E^^^^/E' contains the additional corrections 
for particle showering, causing energy to flow out of or 
into the jet cone. The latter component accounts for 
energy loss from out-of-cone radiation (physics shower- 
ing), leading to a correction of 0.90-1.00 at jet pr >50 
GeV and \y\ <3.0. The equation is iteratively solved us- 
ing Newton's method. The resulting estimate of the jet 
energy is observed to agree with the true E' to better 
than 2% at jet pT > 50 GeV, resulting in less than 0.2% 
uncertainty on the jet response i?had- 

Another issue in using the MPF is related to photon 
identification. To have a clean 7+jet sample in data, only 
CC photons are used with tight selection criteria. How- 
ever, in some jets a large fraction of their transverse mo- 
mentum is carried by photons from tt*', 77, or decays, 
which form a sample of "electromagnetic" jets ("EM- 
jets"). If these photons are sufficiently close together, 
and there is little activity around the photons, the jet can 
mimic an isolated single photon typical for 7-l-jet events. 
Because the cross section for 7-l-jet events is ~3 - 4 orders 
of magnitude lower than that of dijet events, these EM- 
jets contribute a significant background for true 7+jet 
events. An artificial neural network (ANN) is trained 
to discriminate between photon and EM-jets |23j| using 
input variables based on the shape of the calorimeter 
shower and measurements of charged particle tracks in 
the vicinity of the photon candidate. The distribution 
of the photon ANN output for the simulated photon sig- 
nal and for the EM-jet background samples are fitted to 
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FIG. 7: (color online) Extrapolation of the jet energy response 
in the CC at high E' using the reseated MC (see the main 
text) and a fit to the MC points. The dotted line shows a 
simple quadratic logarithmic fit to data for comparison with 
the tuned MC results displayed by the solid line. We also 
display in the bottom panel the relative difference between 
both curves and the statistical uncertainty on the fit to the 
reseated MC in dashed lines. 



the data for each E' and r/ bin using a maximum likeli- 
hood optimization to obtain the fractions of signal events 
in the data. To reduce the uncertainty in the jet en- 
ergy scale due to contamination from background in the 
7-|-jet events, the difference in the response determined 
from real 7-l-jet and dijet events, where one of the jets is 
misidentified as a prompt photon, is estimated using MC 
and applied as a correction based on the estimated purity 
of the selected photons in data. The jet energy response 
after all corrections as a function of E' in the CC is given 
in Fig. [71 The main uncertainty is due to the uncertainty 
on the photon energy scale, which is on the order of 0.5% 
at E' « 20 GeV and 0.8% at E' « 500 GeV. The choice 
of fragmentation model used in Pythia 24] was an ad- 
ditional source of systematic uncertainty on the photon 
purity [111. 

The statistics of the 7-fjets sample limits the direct 
response measurements in the CC to E' < 350 GeV. 
The measured energy response in this region must be 
extrapolated to the highest jet energies at ~ 600 GeV. 
To avoid a statistical uncertainty of more than 2% at 
high-pT in the CC, MC models are used to constrain the 
high-pT response. For this purpose, the measurement of 
the response in 7-l-jet events in the MC is rescaled to 
the measurement in data by modifying the response of 
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jet pt response in the CC: photon energy scale, photon iden- 
tification, fragmentation, and PDF. 



the calorimeter for single pions in MC. Figure [7] shows 
the measured response for jets in data compared to the 
rescaled MC prediction and to a quadratic fit in logE'. 
The uncertainty in the fragmentation model for the high 
E' extrapolation is estimated using the differences be- 
tween the Pythia and Herwig [2^ generators after 
turning off the underlying event modeling. This leads to 
a systematic uncertainty of about 0.8% at E' — 600 GeV. 
The systematic uncertainties related to PDFs (especially 
due to the uncertainty on the gluon fraction in the pro- 
ton) are about 0.2%. 

The total uncertainty on the jet pt response as a func- 
tion of E' is given in Fig. jS) The dominant uncertainty 
comes from the photon energy scale. The uncertainty due 
to photon identification is related to the uncertainty on 
the sample purity and contributes mainly at E' energies 
below 50 GeV. 



2. Pseudorapidity dependent corrections 

The purpose of the 77-dependent corrections is to equal- 
ize the jet response everywhere as a function of pseudo- 
rapidity in the calorimeter after the jets are corrected for 
offset effects. The DO calorimeter is inter-calibrated at 
the cell level as a function of the azimuthal angle 4> by 
equalizing the response of the calorimeter in dedicated (jj- 
symmetric data samples. This yields a jet response that 
is independent of 0, so only the rj dependence of the re- 
sponse needs to be corrected. The 77 dependence of the 
response is mostly due to the changing calorimeter de- 
tector elements, especially in the ICR, different amounts 
of passive material and the varying angle of incidence 
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The sharp features are due to changes in the detector structure, moving from the central to forward calorimeters. 



with jet ry. The 7y-dependent corrections r/) nor- 

malize the response at forward pseudorapidities to that 
measured in the CC (Rcc)- This leads to the definition 

F^{E,ij) = R{E,r^)/Rcc{E), (8) 

where R{E, rf) is the response of the detector for a jet 
of energy E, located at detector pseudorapidity 77. We 
use both dijet and 7+jet samples to determine F^. The 
dijet sample provides high statistics and high reach in jet 
energy for the forward region. One of the jets is required 
to be central and the response measurement is binned 
in terms of the pr of the central jet (using the dijet E' , 
defined as in Eq. [6] where the photon is replaced by the 
central jet) after correcting for the offset and calorimeter 
response. This binning leads to a resolution bias, which 
is later corrected. 

The ?7-dependence of the response, F^, is fitted using 
a quadratic-logarithmic function of E' 

p.w , _ Po(^) +Pi{v) HE') +P2{ri)\v?{E') 



where the pi are fitted as a function of detector rj. The 
Fri and p.^'s are given in Fig. [9l(a-c). As an example 
of the data used in this fit, we give in Fig. [TU] the 77- 
dependent corrections for two bins in rj for the dijet and 
7-1-jet samples. Although the correction factors depend 
on the sample (7-)-jet, dijets), we can remove this depen- 
dency by scaling the dijet correction in the overlap region 
between the CC and the EC by an energy-independent 
factor i?scaic: 

-Rscaic(fy) = 1 + 91 ln[cosh(r7)] + 172 ln^[cosh(?7)], (10) 

where qi and (72 are two parameters fitted to data and 
the result is given in Fig. jHl^d) . This functional form is 
motivated by phenomenological studies of the difference 
in the jet responses measured in 7-l-jet and dijet samples, 
as discussed in the next section. 

The jet pt resolution is worse than the 7 pr resolu- 
tion. Due to the steeply falling inclusive jet cross section, 
more jets migrate into a given px bin from lower px than 
from higher px, giving rise to a px bias compared to the 
particle level. The effect of this resolution bias is taken 
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FIG. 10: Fits of _Fi, in 7+jet and dijet data for two different regions in as a function of E' . The central fit values and the 
uncertainty band are displayed on the figure. 



into account in the final measurement of the jet px re- 
sponse versus rj using the CC jet px resolutions obtained 
from dijet events as described in Sec. IIXI In particu- 
lar, the jet transverse momenta in dijet events in the CC 
are a priori perfectly balanced on average by definition 
[Fj^{E',r] = 0) = i?scaie(?7 = 0) = 1], which provides a 
strong constraint for the bias correction. 

With the application of the dijet-specific scale factor 
and resolution bias corrections we obtain systematic un- 
certainties in the 77-dependent corrections that are less 
than 1% for \tj\ < 2.8 as illustrated in Fig. [Til The 
leading systematic uncertainty is from the average resid- 
uals of the fits for and is estimated to be 0.5% for 
0.4 < \y\ < 2.4 and constant versus energy. This residual 
accounts for the scatter of the data points around the 
central fit and covers possible variation in the shape of 
the fit function. The uncertainty due to the resolution 
bias correction is of the order of 0.5% at jryj = 2.0 and 
reduces to zero at 77 = 0. 



3. Dijet specific response 
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FIG. 11: Relative uncertainties on the r;-dependent correc- 
tions as a function of jet detector rapidity. 



The methods presented so far allow for a precise mea- 
surement of the MPF response in the CC for the 7-l-jet 
sample. However, the response for dijet and 7-l-jet events 
is different. Figure [T^] displays response for the quark 
and gluon initiated jets measured in MC simulations af- 
ter rescaling the single pion response to data. The gluon- 
initiated jets have a lower response than quark-initiated 
jets because they have on average higher particle multi- 
plicity with softer particles. The soft particles lead to a 
lower jet response due to the falling single pion response 
at low energy. Figure [T3] displays the fraction of gluon- 
initiated jets in MC for 7-l-jet and dijet events. The 
7-1-jet jet energy scale cannot be used directly for the 
measurement of the inclusive jet cross section, because 
this sample is strongly dominated by dijets. This effect 
also explains the differences we observe in Fig. [TU] for the 



/^-dependent corrections in 7-l-jet and dijet samples. The 
difference observed in i^,, versus E' at fixed 77 is due to 
the different amounts of quark and gluon jets in the sam- 
ples. The gluon versus quark fractions depend primarily 
on energy (not pt or 77) which leads to a correction factor 
dependent on cosh(?7). Once this difference is taken into 
account, it is possible to combine both samples to fit F^. 

To calculate the relative difference in response between 
7-l-jet and dijet samples in the CC, we first scale the sin- 
gle pion response in MC to reproduce the measured jet 
response in the 7-l-jet data. The measurement from data 
of the absolute jet response in 7-l-jet events in the CC is 
then scaled to its dijet equivalent. The dijet Ty-dependent 
corrections are obtained from a global fit to 7-l-jet and 
dijet data, which accounts for the sample-dependent re- 
sponse. 
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FIG. 12: (color online) Quark- and gluon-initiated jet re- 
sponses and their relative differences for CC jets as a function 
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FIG. 13: (color online) Fraction of gluon initiated jets in 7+jet 
and dijet events in the CC. 



The differences between the dijet response used in this 
analysis and the 7-l-jet response used in most other anal- 
yses are contained in the 7y-dependent scale factor F^i and 
the ratio of tuned MC responses i?dijct/7+jct at 77 = 



RdijctiE' , rj) — -Rdijct/7+jot 



R. 



■7+jot 



(11) 



The ratio between the dijet and 7+jet responses 
^dijct/7+jct is in practice given by the information pre- 



sented in Figs. [12] and [13] and can be expressed using 
the responses for the gluon- and quark-jets^ (^?giuon and 
^quaik) and the fractions of gluon-jets in the dijet and 
7+jet samples (/gi'^°*„ and f^^^n) 



'dijot/7+jct 



f dijet 
luon/„iuon 



-^quark(l 



Y-dijct \\ 
J gluon// 



(-Rgluon/j 



7+jct 
uon 



^quark(l /Jluon)) 



(12) 



C. Showering correction 

Jets are extended objects and deposit their energy over 
a wide area in the calorimeter. When the cone algorithm 
is used, some of this energy is deposited outside the jet 
cone due to interactions with the magnetic field and pas- 
sive material. This is called detector showering and needs 
to be taken into account in the jet energy scale determi- 
nation. In addition, part of the energy of the incident 
parton is lost outside the jet cone because of hadroniza- 
tion and the finite size of the jet cone. This is called 
physics showering and is taken into account in the en- 
ergy scale correction to the particle level. 

The determination of the showering corrections re- 
quires a good understanding of the transverse jet en- 
ergy profile. In a dedicated study, the cell-level in- 
formation from MC is kept to generate energy den- 
sity profiles as a function of the distance Ai? = 
\/ {Uparticle Vjet )^ + {(f'parUcie - 4>jet)'^ between the par- 
ticle and the jet axis for particles originating from inside 
the particle jet, from outside the jet, and from offset due 
in particular to pile up or additional interactions in one 
bunch crossing. The sum of these profiles is fitted to the 
measured energy profile in data to account for possible 
response differences between data and MC. The energy 
profiles are created by summing the energy in the cells at 
a given radius from the cone axis. The profiles are calcu- 
lated for back-to-back 7+jet events and show the jet core 
at Ai? around and the photon contribution at Ai? w tt. 
The energy density in the range i?coiic < Ai? < tt is pri- 
marily offset energy. Figure [14] shows an example of the 
showering profiles in MC without any zero bias event 
overlay (i.e. with only the underlying event and no off- 
set). It gives the average energy in a given rapidity and 
transverse energy bin coming from inside and outside the 
jet as a function of the distance Ai? in rapidity and az- 
imuthal angle from the center of the jet. The MC de- 
scribes the data when both the energies inside and out- 
side the jet are considered. 

The estimate of the showering correction S for 7+jet 
events in MC and data is obtained by comparing the en- 
ergy deposited by all particles inside the calorimeter jet 
cone Y.AR=o Ein+ Z^Afl^o -^out Originating from inside 



^ We note that nearly all the quark-initiated jets come from light 
quarks. 
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FIG. 14: (color online) Jet energy profiles as a function of 
distance from the jet axis AR for MC and data used to com- 
pute showering corrections. The data are corrected for offset 
energy from noise and additional pp collisions and are com- 
pared to MC jets without offset (jet) and contributions from 
the underlying event (not-jet). We note the good agreement 
between data and the sum of energy contributions from inside 
and outside the jet in MC. 



or outside the particle jet to that from the original par- 
ticle jet J2^R=o using the fit-weighted templates 



J2ar=o ^' 



Afl=0 ^out 



EOO p 
A_R=0 ^" 



(13) 



where E-m and -Eout are the energies coming from inside 
and outside the jet. To take into account any potential 
bias in the method, the final value of the showering cor- 
rection in data is computed as 



data 



data 



otruc 

Smc 



(14) 



where the true showering is directly available in 

MC. This bias correction amounts to less than 0.3%. 

While the showering templates are measured in en- 
ergy, the applicable quantity for the cross section mea- 
surement is jet pt- When mapping the showering tem- 
plates to Pt the deposits in rapidity are weighted by 
cosh(j/o)/ cosh(?/i), where yo is the cone axis and yi is 
the rapidity of the energy deposit. As a result of this 
weighting, the effects of showering in px are generally 
suppressed relative to energy showering. This can also 
tilt the jet toward y = and cause a net increase in the 
jet Pt, leading to S'data > 1. The differences between 
energy and pr showering can be up to (1-2)% over the 
kinematic region of the cross section measurements. 

The last step of the showering correction is to make 
the transition from 7-l-jet to dijet events. This remain- 
ing correction is computed directly using the differences 
in showering in 7-l-jet and dijet MC. The final jet px 
showering corrections are given in Fig. 1151 
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FIG. 15: (color online) Jet showering corrections shown as a 
function of jet pr for different regions of jet pseudorapidity. 



The uncertainties on the showering correction are less 
than 1% of the overall correction factor at pT > 50 GeV. 
The main sources of uncertainty come from the difference 
between data and MC in the single pion response at low 
Pt, the quality of the fits of MC templates to data, and 
the description of the underlying event determined by 
varying Pythia tunes for Tevatron data at higher pT- 



D. Potential biases in the method 

1. Topology bias of the MPF method 

The MPF method balances a photon or a central jet 
against the full remaining hadronic recoil, but the mea- 
sured MPF response is interpreted as the response of the 
probe jet. The precision of this interpretation may be 
biased because the hadronic recoil includes particles not 
related to the probe jet, for example, particles coming 
from soft gluon radiation. These additional particles are 
generally softer than those in the core of the jet and are 
expected to lower the response of the recoil with respect 
to that of the jet core. 

In the case of the energy measurement, an additional 
bias is caused by the systematic mismeasurement of the 
jet rapidity, because the MPF method is inherently based 
on balancing pT- As we will see in the following, the ra- 
pidity bias is particularly large in the ICR, where the 
absolute rapidity is systematically underestimated and 
causes a corresponding increase in the MPF response: 
the same calorimeter energy now corresponds to higher 
Pt- Since the bias versus energy has a non-trivial ra- 
pidity dependence and the cross section measurement is 
performed as a function of pT, we derive and apply topol- 



16 



u- 1.05 

Q. 



^ 1.04 



I/) 

o 
o 

Q. 

,o 





R =07- 
f^cone 


- 


-hi = 0.0 : 




--hi = 0.8 : 




hi = 1-6 - 




-.-hi = 2.4 






50 100 


200 300 ,,, 
(GeV) 



o 1 



1.03 



1.02 



1.01 



1.00 



FIG. 16: (color online) Topology bias in the MPF method for 
jet pt response. 



ogy bias corrections as a function of jet px- 

The bias of the MPF response is determined in tuned 
MC by comparing the MPF response to the true response 
defined at the particle level. The result for the px re- 
sponse is shown in Fig. [161 This bias is about 1-2.5% 
for the different rapidity regions with little px depen- 
dence (< 0.5%) for PT > 50 GeV. The MPF response 
bias is quite small since the method is based on the px 
balance and the cone size of R = 0.7 is large enough 
to contain most of the hadronic recoil in the absence of 
additional soft non-reconstructed jets. The bias is signif- 
icantly larger, 2-4%, for R — 0.5 jet cones. The system- 
atic uncertainty on the MPF method bias is computed as 
the difference between the 7-l-jet and dijet samples and 
found to be of the order of 0.1%. 



2. Zero suppression bias 

An off-line zero suppression further suppresses the en- 
ergies of calorimeter cells in order to reduce the amount of 
noise, in particular in the coarse hadronic section, that 
can contribute to jet energies. The algorithm used for 
this zero-suppression retains calorimeter cells if their en- 
ergy exceeds the average baseline noise by 4cr, where a 
represents the measured standard deviation of the noise 
for a given cell. Neighboring cells are also retained if 
their energy exceeds a threshold of 2.5a. 

The zero-suppression algorithm produces a small pos- 
itive noise offset contribution because of the asymmetric 
zero suppression (negative energies are never kept). For 
cells with high enough real energy deposits, as within 
the jet core, the zero-suppression produces no net offset, 
and positive and negative noise offset contributions are 
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FIG. 17: Bias in the measurement of jet energy, ko, due to 
zero suppression effects on the offset correction, shown as a 
function of p'x for central jets and different number of recon- 
structed vertices. 



expected to cancel. Conversely, the energies measured 
for particles incident on the calorimeter, including those 
from uranium decay, are reduced by the zero-suppression 
when cells are below threshold. Therefore the average 
offset within a jet is different from the offset outside of a 
jet which we measure using zero bias and minimum bias 
events. 

The true offset is increased inside the jet environment 
compared to the average energy density measured outside 
jets in zero bias and minimum bias events. The correction 
factor for the true offset inside a jet, fco, is defined as 



kn = 



(EtLsi^^o ZB overlay)) 
EiZ^ZB overlay) - O 

mca 



(15) 



jjjpjjs is the energy of a reconstructed jet and 



where 

Omcas is the measured offset correction described in 
Sec. IV Al The same MC events are reconstructed with 
and without zero bias event overlay (offset). The zero 
bias event sample was collected without any calorimeter 
zero suppression so that its effect can be studied in de- 
tail. Figure [17] shows the effect of zero-suppression on 
the offset correction for jets in the CC. For jet p'x > 50 
GeV, where p'x = i?'/cosh rjjct, the resulting bias on jet 
energy varies between 5% at low pip and 2% at higher pip. 

The bias in offset is almost fully canceled by an oppo- 
site bias in the MPF response, defined as 



kn = 



(^had(no ZB overlay)) 
(i?had(ZB overlay)) 



(16) 



because the increased offset inside the jet increases the 
in the photon direction. This artificially increases 
the estimated MPF response (see Eq. 5). Only the ratio 
ko I kfi is therefore relevant for the final bias correction 
due to the zero suppression bias. The combined bias is 
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found to be less than 0.5% for jet pr > 50 GeV in all 
rapidity bins, largely cancelling the topological bias, and 
approaches zero at high pT, as shown in Fig. [THl 
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FIG. 19: (color online) Rapidity bias obtained in MC for 
different jet energies E. The curves cover the range pr > 
30 GeV and lyjct] < 2.4. 



3. Rapidity bias 

Since the inclusive jet cross section is measured in bins 
of rapidity, we checked for any potential bias in the recon- 
struction of jet rapidity using the simulation, as shown 
in Fig. [ini The rapidity is generally biased towards the 
central calorimeter, with the largest deviations observed 
in the ICR. This is attributed to detector effects in the 
ICR in addition to the jet cone algorithm itself. The 
absolute effect on the inclusive jet measurement is small 
compared to the effect of jet px calibration. 



E. Final jet energy scale corrections and 
uncertainties 

Figure [20] shows the jet energy scale corrections as a 
function of jet pr for central and forward rapidity, and 
as a function of jet rapidity at low and high jet pr- The 
corrections range between 1.2 and 1.8 for the kinematic 
range of the cross section measurement. The response 
correction is by far the largest one, while the showering 
correction starts to be noticeable at large rapidity. At 
high rapidity, the actual angular distance for each Ar] 
bin is small, while the radius of the showering is slightly 
increasing due to the increasing energy of the jet at fixed 
Pt as one goes forward. The total correction is computed 
using Eq. 3. The combined effects of the uncertainties as- 
sociated with each component of the correction are sum- 
marized in Fig. [5T]as a function of jet pt for central and 
forward rapidity, and as a function of jet rapidity for low 
and high jet pT — high jet energy corresponds to low pT 



at high rapidity. 

The corrections do not show a significant dependence 
as a function of jet rapidity except in the region of the 
ICR. The uncertainties vary between 1.2 and 2.5% for the 
kinematic range of the cross section measurement and are 
dominated by the uncertainties of the jet response. The 
uncertainties obtained in the CC and for jet pT ~ 100- 
500 GeV are the smallest ones obtained by any experi- 
ment operating at a hadron collider. These uncertainties 
do not depend strongly on jet pseudorapidity and pt- 

F. Closure tests 

The aim of the closure tests is to verify the accuracy of 
the jet energy scale correction using either MC or data 
and to evaluate the remaining difference as an additional 
systematic uncertainty related to the method. As an ex- 
ample, one test is to use the full method in MC and 
to compare the results with the particle level jet energy. 
The direct closure variable is defined as 

/ TT^corrX 

\ , (17) 

/ T-iparticle \ 

where Ej°" is the corrected jet energy and £'P^'^''"=^<= jg 
the energy of the closest particle jet matching the re- 
constructed jet within AR < i?cone/2. Results from the 
direct closure test are shown in Fig. [22] in two regions 
of jet rapidity. We note that we obtain consistency of 
the method within statistical uncertainties (D is close to 
unity within less than 1%) and no additional systematic 
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FIG. 20: (color online) Jet energy scale corrections as a function of jet pr for (a) central and (b) forward rapidity, and as a 
function of jet rapidity for (c) low and (d) high jet pt- 



uncertainty is introduced. Closure tests using data are 
performed relative to MC by comparing ratios of fully 
corrected jet energies < EXtT"^""^ > / < ElfJ'"=*'"^ > in 
fixed regions of E' and ry. Again we find good agreement 
within the expected uncertainties of the jet energy scale. 



VI. TRIGGERING ON JETS 

In this section, we briefly describe how we determine 
the absolute jet trigger efficiencies. Two different samples 
based on jet or muon triggers are used. The DO trigger 
system is composed of three consecutive levels called LI, 
L2, and L3. At LI, a single jet trigger typically requires n 
calorimeter trigger towers above a given threshold, where 
a trigger tower is defined by the hardware summation of 
energies in 2 x 2 calorimeter towers. The trigger towers 
are read out separately from the precision calorimeter 
electronics via a fast digitizer and are used in both LI 
and L2 triggers. All events used in this analysis are re- 
quired to pass a trigger designed to fire if a single jet 
with pt > 50 GeV is in the event. For instance, the 
65 GeV single jet trigger requires the presence of three 
calorimeter towers with a transverse momentum above 5 
GeV. This requirement is often satisfied by the presence 



of trigger towers belonging to different jets, ensuring high 
trigger efficiency. In most of these events, there are two 
high-pT jets in the event or more than two low-px jets, 
which ensures that the event passes the LI threshold. A 
detailed analysis shows that the LI single jet efficiency is 
more than 98% for the full kinematic range of our mea- 
surement, which is corrected for the residual inefficiency. 
At L2 we perform a clustering of the trigger tower en- 
ergies and apply a threshold based on the px of highest 
energy cluster. Seven L3 triggers corresponding to un- 
corrected L3 jet Pt thresholds of 8, 15, 25, 45, 65, 95, 
and 125 GeV are used in the analysis. The highest-pr 
L3 trigger was never prescaled during data collection. In 
Fig. [531 we show the jet cross section before any unfold- 
ing corrections as a function of jet pt for the different 
jet triggers for two domains in jet rapidity \y\ < 0.4, and 
2.0 < \y\ < 2.4. 

The first method used for computing the jet efficiency 
is to obtain the relative jet trigger efficiency with respect 
to the lower pT jet trigger. For instance, the ratio of the 
95 and 65 GeV triggers is shown in Fig. [Ml For this 
purpose, we plot the ratio of the number of events that 
pass the 95 GeV trigger to those that pass the 65 GeV 
threshold as a function of jet pt after jet energy scale 
and vertex efficiency corrections to cancel known lumi- 
nosity dependencies as discussed in Sec. VII. The ratio 
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FIG. 21: (color online) Jet energy scale uncertainties as a function of jet pr for (a) central and (b) forward rapidity, and as a 
function of jet pseudorapidity for (c) low and (d) high jet pr- 
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FIG. 22: Closure test of jet energy scale in MC for (a) \ridet\ < 0.4 and (b) 0.4 < \ridet\ < 0.8. The band outlined by the solid 
curves corresponds to the uncertainties in the extraction of jet energy scale in MC which are mainly statistical. 



is scaled by the relative integrated luminosities of these 
triggers to account for the different prescales. When this 
ratio reaches 1, the 95 GeV threshold trigger is 100% ef- 
ficient with respect to the 65 GeV one. A fit to this ratio 
gives the different thresholds for which the triggers are 
fully efficient (>99.9%). The jet energy scale corrected 
Pt at which each trigger becomes fully efficient is given in 
Table |T| These thresholds take into account the px bin- 
ning used in the analysis and can be significantly higher 



than the minimum usable threshold. We note that this 
method does not allow us to obtain the absolute trigger 
efficiency since it gives all efficiencies with respect to the 
lowest 8 GeV pr trigger as a reference. 

A second method is used to measure the absolute sin- 
gle jet trigger efficiency. It uses independent muon and 
minimum bias triggers. The minimum bias trigger only 
requires energy deposits in the luminosity monitors. As 
its name indicates, it shows very little selection bias and 
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is ideal for trigger studies. Unfortunately, the sample 
collected during all of Run II at the Tevatron at 0.5 Hz 
only yields statistics adequate to study jets below 70 GeV 
using this trigger, and this method does not allow explo- 
ration of the high px jet trigger efficiency. For this reason, 
inclusive muon triggers without any calorimeter require- 
ments are also used. This allows us to check what fraction 
of the offline reconstructed jets in muon triggered events 
pass the calorimeter jet trigger requirement, providing a 
direct estimate of trigger efficiencies up to 400 GeV in the 
CC. The conclusion of this study is that all jet triggers 
are more than 98% efficient above the thresholds defined 
above, and the residual inefficiency is determined to a 
precision of better than 1%. Both methods to obtain the 
trigger efficiencies are useful since the muon triggers have 
a tendency to enrich the inclusive jet samples in b and 
c-jets where the b and c quarks decay leptonically, which 
might lead to different trigger efficiencies as a function of 
jet pT- 
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TABLE I: Jet energy scale corrected pr in GeV at which each 
L3 trigger becomes fully efficient in different jet y bins. 



VII. EVENT SELECTIONS AND EFFICIENCIES 

In this section, we discuss the selections that are used 
to remove background events in the sample. The selec- 
tions fall into three different categories. The event qual- 
ity flags remove events suffering from diverse calorimeter 
noise issues. The vertex requirement selects events with 
a high quality vertex close to the center of the calorime- 
ter to improve the jet pr and y measurements and to 
reduce the background from cosmic ray events. The J^j^ 
requirement is designed to remove the remaining cosmic 
ray background, especially at high jet px- 



A. Event quality flags 

Event quality flags ensure that the subdetectors used 
in the analysis were working properly when the data were 
collected. Calorimeter event quality flags allow removal 
of events showing coherent pedestal shifts in the analog- 
to-digital converters, parts of the calorimeter not cor- 
rectly read out, or high coherent noise. This is espe- 
cially important for high-py jets which can originate ar- 
tificially from noisy towers in the calorimeter. Note that 
the vertex and l^x requirements also remove most of these 
events. The inefficiency induced by the calorimeter event 
quality flag rejection is estimated using an independent 
sample whose trigger is known to be unaffected by the 
calorimeter problems, the zero bias trigger. The ineffi- 
ciency is calculated to be (3.2 ± 1.0)% where the 1.0% 
uncertainty covers the time and luminosity dependence 
of the inefficiency. 



21 




100 150 200 250 300 350 400 



(GeV) 



FIG. 24: (color online) Ratio of two consecutive jet triggers 
used to obtain the relative jet trigger efficiency for the 95 
and 65 GeV single jet triggers. The fit of the turn-on curve 
determines the 95 GeV jet pr 99%-threshold to be 130 GeV. 
A higher threshold of 160 GeV that is consistent with an 
efficiency e = 1.00 is used in the final analysis. 



B. Reconstructed vertex requirement 

The vertex selection is based on three different require- 
ments: there must be at least one reconstructed vertex, 
the z-position along the beam line of the primary recon- 
structed vertex must be within 50 cm of the detector cen- 
ter (IzvortcxI < 50 cm), and the number of tracks fitted to 
the vertex has to be at least three to ensure an accurate 
measurement. The z-vertex position requirement ensures 
that the vertex is in the high efficiency tracking region. 
The third requirement rejects vertices originating from 
fake high px tracks. To each reconstructed vertex is as- 
signed a probability that it comes from a minimum bias 
interaction based on the hi{pT) distributions of the tracks 
with pt > 0.5 GeV pointing to the vertex. The vertex 
with the lowest minimum bias probability is selected as 
the primary vertex. 

The efficiency of reconstructing a vertex with at least 
three tracks pointing to it (without the requirement on 
the z-vertex position) is found to be (99.6 ± 0.4)%, inde- 
pendent of jet Pt and y. The observed 0.4% inefficiency 
is consistent with about 0.6% of the primary vertices 
not being reconstructed because of tracking inefficiencies, 
and 0.2% being replaced by a minimum bias vertex. 

The leading inefficiency comes from the requirement 
on the vertex position along the z-axis. The fraction of 
events rejected by this requirement is of the order of 7%. 
To determine the efficiency of this requirement we take 



into account the shape of the luminous region. The lon- 
gitudinal shape of the luminous region is approximated 
by the expression 



(18) 



where the overlap of the proton and antiproton beam 
bunches having Np and Np particles is described with a 
Gaussian distribution of width cr^ in the z direction, with 
a possible offset zqz relative to the nominal interaction 
point. (Txiz) and <Jy{z) represent the transverse size of 
the beam spot and vary as a function of z: 
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(19) 



Here T is either x or y, zqt is the minimum of the /3 func- 
tion describing the beam dimensions near the interaction 
point in direction T and any offset in the x and y direc- 
tions with respect to the nominal interaction point, 7 is 
the Lorentz factor of the beam particles. The emittance 
ct and beta parameter describe the beam dimensions 



at the interaction point, 
integrated to yield 



The parameterization can be 
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where the limits of integration in the denominator come 
from the requirements used in the luminosity determi- 
nation. This parameterization is fitted to minimum bias 
data in the high tracking efficiency region (|zvertex| < 40- 
60 cm) in bins of instantaneous luminosity for several run 
ranges (the changes in beam optics as a function of time 
affect the beam shape as described by the /3* parame- 
ter). The changes as a function of instantaneous lumi- 
nosity are primarily due to the variations in the beam 
parameters during a store. The vertex efficiency varies 
by up to 6% as a function of instantaneous luminosity 
and by up to 4% as a function of the period of data tak- 
ing for a fixed value of luminosity. The parameterizations 
have been determined as a function of time and instanta- 
neous luminosity, and are applied as such on a per-event 
basis. Figure [21] shows the mean vertex efficiency as a 
function of instantaneous luminosity, with the range of 
efficiencies overlaid. The uncertainty on the vertex ac- 
ceptance is estimated to be 0.5% by comparing results 
from fits to minimum bias data at |zvcrtcx| < 60 cm and 
l-^vcrtcxl < 40 cm. In addition, an increased uncertainty 
of 0.4% added in quadrature at high \y\ is introduced to 
account for the possibility of a lower vertex reconstruc- 
tion efficiency. 



C. Missing transverse energy requirement 

A requirement on the missing transverse energy in an 
event is applied to remove the remaining background 
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FIG. 25; (color online) Vertex acceptance for the require- 
ment on the 2- vertex position |zvcrtcx| < 50 cm as a function 
of instantaneous luminosity. The shaded band indicates the 
variation for different running periods. 



from cosmic rays that induce showers in the calorime- 
ter. The cross section for these cosmic ray interactions 
falls much less steeply versus px than the inclusive jet 
cross section, and is typically comparable at pr ~ 400 
GeV. The issue of background from cosmic rays is thus 
more important for high-p^ jets. Fortunately, cosmic ray 
showers deposit most of their energy on one side of the 
calorimeter, have no reconstructed vertex, and produce 
high uncorrected l^rp that peaks at j3T,ioad/^T ~ 1 where 
PT,iead IS the Uncorrected pt of the leading jet of the 
event. These events are fully and efficiently removed by 
requiring pr, lead /I^t > 1-4, when pT.icad < 100 GeV and 
PT,icad/^T > 2.0, when pT.iead > 100 GcV. Figure M 
shows the distribution of pt, \oa.d / for the high-p^ jet 
trigger with px > 200 GeV, with the selected events at 
PT.\cad/$T > 2.0 shown by the shaded region. A spike 
coming from cosmic ray events is visible at 1. An upper 
limit of 0.4% is estimated on the inefficiency of the l^x re- 
quirement and used as an uncertainty, but no correction 
is applied. This upper limit is based on studies of fits of 
distributions like the one in Fig. 1261 and track-matching 
inefficiency for jets since cosmic ray events are usually 
out-of-time with the tracking read-out. 



VIII. JET IDENTIFICATION REQUIREMENTS 
AND EFFICIENCIES 



The jet identification requirements are designed to re- 
move instrumental backgrounds such as jets formed from 
sources of transient noise in the calorimeter and also 
physics background from electrons and photons. The jet 
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FIG. 26: (color online) Distribution of PT.icad/^T for J^t 
events with leading jet pt > 200 GeV. A peak from cosmic 
ray background is visible around 1. The shaded region shows 
jets passing the PT,ioad/^T requirement. 



requirements are based on the fractions of jet energy de- 
posited in the electromagnetic calorimeter (EMF) and in 
the coarse hadronic calorimeter (CHF). EMF< 0.95 is re- 
quired to remove overlaps between jets and electromag- 
netic objects, i.e. electrons and photons. This retains 
true jets with a 99% efficiency. A lower limit on EMF 
(either or varying between 0.03 and 0.05 depending on 
the pseudorapidity region in the calorimeter) as well as 
an upper limit on CHF (varying between 0.4 and 0.6) re- 
moves jets that are formed predominantly out of noise in 
the hadronic calorimeter. An additional requirement, LI 
confirmation, is based on the ratio of the px as measured 
by the LI trigger system and as measured by the preci- 
sion read-out. It is required to be above 0.5 for jet pt < 
80 GeV, and there is no requirement for higher px jets. 
This removes jets formed out of noise, for example due 
to coherent noise in the precision readout electronics. 

The jet identification efficiencies are determined using 
a data driven method. This method uses track jets which 
are jets built with a cone algorithm using charged particle 
tracks instead of calorimeter energy clusters. We select a 
leading px tagged object, which in this case is a photon 
or a track jet associated with a good calorimeter jet, and 
a probe object, which is the leading track jet that is back- 
to-back in (f) with the tag object. Events with additional 
track jets are vetoed to ensure that the leading objects 
are balanced in px- The reconstruction efficiency is de- 
fined as the fraction of probe objects with a calorimeter 
jet found within the 0.7 jet cone, and the jet identification 
efhciency is the fraction of those calorimeter jets pass- 



23 



O 

o 1 

0) 

§ 0.98 

CD 



0.96 
0.94 



^ D0, 0.70 fb" 



i T T T 1 



ideti 



< 0.4 

Dijet tag-and-probe 
Central fit 

Systematic uncertainty 



50 



100 



150 



200 



250 



300 



(GeV) 



FIG. 27: (color online) Jet identification efficiencies obtained 
for the dijet sample. Dashed lines indicate the systematic 
uncertainty. 



ing the jet identification requirements. The data driven 
method has been used for three different samples: dijet, 
7+jet and, Z+^et, which all lead to the same result. The 
central value for the jet identification efficiency shown in 
Fig. [27] is taken from the dijet sample. The efficiency for 
Pt > 50 GeV, where we perform the measurement of the 
inclusive jet px cross section, is 99% in all calorimeter 
regions except in the region 0.8 < \y\ < 1.2 where it is 
about 98%. 

Because the data driven method is used for calorime- 
ter jets that are independently identified as track jets, 
we also directly measure the efficiencies by computing 
the fraction of events removed by each jet identification 
requirement individually after applying all other require- 
ments in the inclusive jet sample. This method assumes 
that each jet identification cut removes only good jets. 
The efficiencies described above are found to be in good 
agreement with those from the tag-and-probe method. 



A. Dijet asymmetry 

The jet px resolutions are determined starting from 
the dijet asymmetry 



Pt.i -Pt,2 
Pt,i + PT,2 



(21) 



computed in a pure dijet sample with no additional jet 
identified, where pt,i and pT^ are the pt of the leading 
and second-leading jets and the two leading jets are ran- 
domly assigned an index of 1 or 2. Both jets are required 
to be back-to-back with > 3.0 to avoid any large ef- 
fects from QCD radiation. The RMS of the asymmetry 
distribution is directly proportional to the jet pt resolu- 
tion 



1 (Jprj, 

V2 PT 



(22) 



if the jets are in the same y region to ensure that the pt 
resolution of both jets is the same. To characterize the 
PT dependence of the resolution for a single jet, ga is 
measured in bins of pt = {pt,i +Pt,2)/2. This method 
can be used directly to measure the jet pr resolution in 
the central region where the statistics are high. However, 
in the forward region, the statistics for forward-forward 
jet pairs is small compared to central-forward jet pairs. 
If one of the jets is in the central region and the other 
in the forward region, it is possible to infer the jet pT 
resolution (Tp^ in the forward region once the resolution 
for jets in the central reference region (Ticf is known 



(23) 



The central reference region used in this study is |?/rof | < 
0.8, with the probe jet binning following the same 0.4 bin- 
ning in rapidity as the rest of the analysis. The asymme- 
try distribution in the central region is shown in Fig. [28] 
for 80 < Pt < 100 GeV as an example and other p-!^* 
bins also show similarly small non-gaussian tails. 




IX. JET TRANSVERSE MOMENTUM 
RESOLUTION 

In this section, we discuss the determination of the jet 
Pt resolution, which is needed for the unfolding of the 
inclusive jet pt cross section. The jet pT resolution is 
determined from data using the dijet asymmetry distri- 
bution, which can be obtained with minimal input from 
MC. This method requires corrections for the presence of 
additional unreconstructed jets (soft radiation), momen- 
tum imbalance at the particle level, and asymmetry bias 
due to non-Gaussian tails. We describe each correction 
needed to obtain the jet pt resolution. 



B. Corrections to the resolution 

The jet pT resolution determined from the dijet asym- 
metry can be affected by physics and instrumental ef- 
fects. The final parameterization of the resolution used 
in this analysis includes corrections to remove biases in 
the measurement as described below. 



1. Soft radiation corrections 

The asymmetry method to compute the jet px resolu- 
tion is biased by the presence of non-reconstructed jets 
in the sample. The pt threshold to reconstruct a jet is 6 
GeV, and requesting the presence of only two jets in the 
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FIG. 28: (color online) Asymmetry distribution for jets in 
the central region with 80 < pr < 100 GeV. The probe jet 
is at \y\ < 0.4, the reference jet at |j/ref| < 0.8. The two 
lines display the result of a Gaussian fit and a Gaussian with 
smeared exponential tails (see Sec. lIXC|) . 



sample to compute the asymmetry does not ensure the 
absence of jets with pr below 6 GeV. The corrections 
for such soft radiation are determined directly in data. 
We compute the asymmetry and the jet px resolution for 
different px thresholds for jet reconstruction, namely 7, 
8, 10, 12, 15, 20, and 40 GeV. The jet pr resolution as 
a function of the jet reconstruction threshold is shown in 
Fig.[29](a) for one bin in jet pt and \y\. A linear fit allows 
for extrapolating the jet pt resolution to a threshold pt 
of 0. The soft radiation factor, 



soft 



iPT%t ^ 0)/PT 



7pT (Pt "oft 



6 GeV)/pT' 



(24) 



is studied as a function of the average jet px in each \y\ 
bin as illustrated in Fig. [51] (b) . To better describe the 
low Pt region and limit the statistical fluctuations, the 
dependency of ifgoft versus px is fitted with 

^soft(PT) = 1 - exp(-po -PiPt), (25) 
where po and pi are two parameters of the fit. 

2. Particle imbalance and combined corrections 

The remaining correction needed to obtain the final jet 
px resolution is the particle imbalance correction. Even 
in the ideal situation of only two particle jets and no 
soft radiation, the two jets are not necessarily perfectly 
balanced. In particular, fragmentation effects cause some 
energy and px to be found outside the jet cone. This 
effect is purely related to QCD and is determined using a 
MC simulation. The particle level imbalance is corrected 



for soft radiation using the same method as introduced 
for data 
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where cr^^^ is the resolution evaluated at the particle 

level in the MC and Px'ptfi"^'^ is the px threshold of jet 
reconstruction at the particle level. 

The corrected particle level imbalance (Tmc is sub- 
tracted in quadrature from the soft-radiation corrected 
resolution computed in data (see previous section). 



'MC 



(28) 



The relative correction due to particle level imbalance 
is about (7-9)% in the CC, (2-6)% in the ICR and the 
EC, for Px > 50 GeV. The systematic uncertainties on 
particle imbalance corrections are mainly due to the dif- 
ferences between the Gaussian one standard deviation 
and the RMS of the particle level imbalance distribu- 
tion due to non-Gaussian tails. The RMS is used for the 
central correction. The main non-Gaussian tails in parti- 
cle level imbalance corrections are caused by muons and 
neutrinos, which are not included in the definition of DO 
particle jets. 



C. Final jet px resolutions 

Using the asymmetry method and the various correc- 
tions discussed above, we obtain the jet px resolutions 
shown in Fig. 1301 The measured resolutions are fitted 
with the parameterization 




(29) 



where N is the noise term, S the stochastic term, and 
C the constant term. The values of the parameters are 
given in Table |TT1 These resolutions are used to obtain 
the inclusive jet px cross section as described in the next 
section. 

We note that the resolution is not Gaussian at high px 
even in the central region because of calorimeter punch- 
through: jets at very high px are not always fully con- 
tained in the calorimeter and can deposit energy into the 
muon system. In Fig. [3T] (a), we show the distribution 
of (pT,jot/PT,ptci — 1) — the ratio of the reconstructed to 
the particle level jet px — obtained from MC simulation 
of the detector in the central region of the calorimeter 
at high Px- The ICR also exhibits non-Gaussian tails as 
shown in Fig. [31] (b), which are explained by the chang- 
ing structure of the calorimeter in this region. The non- 
Gaussian tails are modeled using a smeared exponential 
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FIG. 29: (a) Jet pr resolution extrapolated to a jet pr reconstruction threshold of GeV (in the 80 < pr < 100 GeV and 
\y\ < 0.4 bin), (b) Soft radiation correction factor as a function of the average jet pr for the 0.4 < \y\ < 0.8 bin. The solid 
lines are the fit result and the dashed lines are the result of the extrapolations. 
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with /i, (T, P, and A as free parameters. The fitted shape 
from MC is scaled by varying the parameter a such that 
folding the distributions for the leading jets with the ex- 
ponential pt spectrum from data results in precisely the 
same RMS of the jet px resolution as observed in data. 
This method can account for any shaping of the non- 
Gaussian tails that takes place due to bin-to-bin migra- 
tions in data. The full MC shape with tuned a is later 
used in the unfolding of the data. 

The uncertainties on jet pT resolution are given in 
Fig.[32]for two bins in rapidity as an example. The uncer- 
tainties come primarily from the statistical uncertainties 
in the fits. An additional component is added to cover 
non-statistical variations between the fit model and the 
data. The total uncertainty coming from the jet pt res- 
olution is (5-10)% over the full kinematic range covered 
by the inclusive jet cross section measurement {pr > 50 
GeV). The leading systematic uncertainty in the central 
region is (4-5)% due to the uncertainties on the particle 
level imbalance corrections. In the ICR, an important 
systematic is due to the uncertainty on the tails in the 
resolution for this region. This systematic is estimated 
by varying the size of the tails by a factor of two, and 
is not included in the RMS of the resolution, but rather 
the resulting variation in shape was used in the unfolding 
procedure for data. Another important source of uncer- 
tainty is taken from the following MC closure test: the 
full resolution measurement using the asymmetry is re- 
done using a full simulation of the DO detector, and the 
difference between the MC input true resolution and the 



result of the method is taken as a systematic uncertainty. 
This amounts to up to about 10% uncertainty in the res- 
olution at Pt — 50 GeV in the forward region. 
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TABLE II: Parameters of the fits to the jet pr resolution 
versus pr for data. The noise term is fixed to the MC value 
with an uncertainty of 1 GeV since it is not constrained by 
the data. 



X. JET RAPIDITY RESOLUTION 



Compared to the jet pt resolution, the rapidity resolu- 
tion is a small effect which is determined using a MC sim- 
ulation of the detector. The bin width in y is much larger 
than the y resolution and bin-to-bin migrations only oc- 
cur at the bin edges. To unfold the effect of the rapidity 
resolution, a smooth parameterization of the resolution 
as a function of y is used. The result of the parameteri- 
zation of the y resolution in different px bins is shown in 
Fig. El 



26 



.0.20r 



D0, 0.70 fb"^ 

|y| < 0.4 

Rcone = 0.7 

X2/NDF = 2.3 / 11 




30 40 



100 



200 300 



.0.20r 



(GeV) 



D0, 0.70 fb"^ 
0.8< |y| < 1.2 




.0.20r 



^H0.18: 



0.16: 
0.14: 
0.12: 
0.10: 
0.08: 
0.06: 

.0.20r 




D0, 0.70 fb"^ 
0.4< |y| <0.8 

Rcone = 0.7 

X2/NDF = 5.0/11 



(b) 



30 40 



100 



200 300 



(GeV) 



0.18 

O.I67 |\I 
0.14 

0.12; 
0.10; 

0.08 



D0, 0.70 fb"^ 
1.2< |y| < 1.6 




(d) 



0.06 



X2/NDF = 10.7/10 



30 40 



100 



200 300 



Q. 
D 



.0.20r 
■0.187 
0.167 
0.147 
0.127 
0.107 

0.08- 
0.06- 



P^ (GeV) 



30 40 



100 



200 300 



D0, 0.70 fb"^ 
1.6< |y| <2.0 

Rcone = 0.7 

X2/NDF= 1.2/8 



(e) 




.0.20r 



P^ (GeV) 



Q. 
D 



D0, 0.70 fb"^ 
2.0< |y| <2.4 

Rcone = 0.7 

x2/NDF = 8.2/8 




30 40 



100 



200 300 



P^ (GeV) 



30 40 



100 



200 300 



P^ (GeV) 



FIG. 30: (color online) Jet pr resolution determined in data for the six rapidity regions. The solid curves are the results of the 
fit. The fit uncertainty is given by the shaded band. 
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FIG. 31: (color online) (a) Example of punch through for < 0.4 at high jet pr for 300 < pr < 400 GeV and \y\ < 0.4. (b) 
Example of tails of the jet pr resolution in the ICR for 160 < pr < 200 GeV and 0.8 < jyj < 1.2. The two curves are the result 
of the Gaussian fit and of a Gaussian plus exponential tails. 
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(color online) Relative statistical and systematic uncertainties on jet pr resolution for (a) \y\ < 0.4 and (b) 2.0 < 
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FIG. 33: (color online) Rapidity resolution (RMS) as a func- 
tion of y in different jet pr regions. 



XI. UNFOLDING 

In this section, we describe the method used to unfold 
the data as a function of jet px and y. As we already 
mentioned, the main smearing effect is due to the jet pr 
resolution while the y smearing is only a second order 
effect. The steeply falling jet pt cross section convoluted 
with the jet px resolution leads to an increase of the ob- 
served cross section as a function of the measured jet pt- 
To unfold the data, we use the so-called ansatz method. 
We start with a functional form for the cross section that 
has only a few parameters, smear it with the jet pr and 
y resolutions, and fit the parameters so that it describes 
the raw cross section measurement before unfolding. 

The ansatz used in each rapidity bin contains a px 
dependence term and an additional rapidity dependence 



PT 



100 Gev) 
exp {-jpt) ■ 



2pT cosh(|y„ 



/9 



(31) 



Here ^/s — 1960 GeV is the center-of-mass energy and 
I Vmin I is the low edge of the bin in absolute rapidity. The 
ansatz is based on phenomenological fits and motivated 
by the parton model [l^]. The exponential term repre- 
sents hydrodynamic production by freezing out particles 
from the quark and gluon sea. The value of 7 is ex- 
pected to be of the order of 0.3-0.6 GeV^^, typical of 
the proton size. The first power term characterized by 
a represents the scaling violations associated with hard 
production. Typical values of a are 4-6 for single particle 
production. The second power term characterized by /3 
represents the kinematic suppression effect at the edges 
of the phase space of particle production. 

The ratios between the data and the smeared ansatz 



are shown in Fig.[34l where the ansatz correctly describes 
the data in all y bins. The unsmearing corrections for 
the PT resolution effects are shown by the dashed lines in 
Fig. [3SJ The unfolding corrections are (10-40)% in the 
CC, (20-80)% in the ICR where the jet pr resolution is 
worse, and (15-80)% in the EC where the jet cross sec- 
tion falls steeply. The highest pT bin (where the unfold- 
ing corrections are the largest) where the cross section is 
measured is chosen so that the cross section measurement 
is still meaningful; the number of events should still be 
sufficient to give a lower limit on the measured cross sec- 
tion at the 95% C.L. (iVtheory/ V A^smeared theory > 1.645). 

Although in some bins most of the events migrate from 
lower Pt, the migrations are well understood and result 
in a relatively small uncertainty compared to the uncer- 
tainty from the jet energy scale. The ansatz unfolding is 
found to be in good agreement with the results using the 
Pythia MC where the cross section is rescaled to data 
and the jets at particle level are smeared according to the 
Pt resolutions obtained in Sec. IXII 

The same ansatz unfolding method can be used to un- 
fold the cross section for effects of the resolution for re- 
solving rapidity, assuming the pt and y resolutions are 
uncorrelated. Since the y resolution is much better than 
the pt resolution, the effects of the y resolution are a 
small perturbation on top of the pT smearing. The fits 
to the unfolded pt spectra (unfolded for pT resolution ef- 
fects only) in neighboring rapidity bins are interpolated 
with respect to rapidity to produce a smooth, continuous 
two dimensional spectrum in pt and y. A final unfold- 
ing is performed to correct for events that migrate into 
neighboring rapidity regions due to effects of the y reso- 
lution. The results of the y unfolding as a function of jet 
Pt in the different y bins are given in Fig. 1351 together 
with the results of the global unfolding corrections in jet 
Pt and y. As expected, the effects of y unfolding are very 
small with respect to the effects of the pt unfolding. 



XII. INCLUSIVE JET px CROSS SECTION 
MEASUREMENT 



In this section, we describe the final result on the in- 
clusive jet Pt cross section measurement applying the 
corrections defined in the previous sections: jet energy 
scale, efficiencies, and unfolding, used in order to com- 
pute the true number of events observed in each pT and 
y bin. The cross section results are given in Fig. [551 in 
the six y bins as a function of jet pT- The data points 
are plotted according to the prescription described in |23| 
and the tabulated data are available from Ref. (29j . 

The method used to extract the cross section is re- 
peated and cross checked using a MC simulation of 
the detector. Events are generated using PYTHiA and 
weighted to match the NLO prediction calculated us- 
ing the CTEQ6.5M PDFs and including nonperturbative 
corrections. The MC events are treated in the same way 
as data, all corrections are rederived using MC events. 
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FIG. 34: Data divided by the ansatz fit witli models for pr and y smearing in tiie six rapidity regions. 
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FIG. 35: (color online) Unfolding corrections in the six rapidity regions as a function of jet pt- The corrections are given for 
Pt, y unfolding separately and combined. 
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FIG. 36: (color online) Inclusive jet cross section measure- 
ments as a function of jet pr in six \y\ bins. The data points 
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and the derived cross section is compared to the input 
cross section to perform a closure test of the measure- 
ment. The results given in Fig. [37] show that the method 
used to extract the cross section works well within the 
statistical uncertainties of the fits to the jet response, jet 
Pt resolution and pt spectrum. These MC uncertainties 
are significantly smaller than the systematic uncertain- 
ties present in data. 

In Fig.|36]the measurement is compared to the predic- 
tion of NLO QCD using the CTEQ6.5M PDF parame- 
terization computed using the nlojet-| — h program and 
FastNLO. The central CTEQ6.5M prediction uses the 
factorization and renormalization scales iip = fxn — pt- 
The alternative scale choices fip = fJ-R — O-Sipr and 
fJ-F = fJ-R = 2pT are used to estimate the theoretical 
uncertainty on the higher order corrections. 

The NLO PQCD prediction is corrected for non- 
perturbative effects to connect the parton level jets pre- 
dicted by theory to the measured particle level jets. 
The leading non-perturbative corrections are hadroniza- 
tion and underlying event which partially cancel. An- 
other small correction is the exclusion of muons and 
neutrinos from the definition of the particle jets. The 
muon/neutrino energy loss is not corrected by the JES 
procedure using the MPF method in 7-l-jet events. The 
MC corrections have been obtained using Pythia v6.412 
with parameters for tune QW [s^l obtained by tuning 
Pythia to reproduce CDF Run II data. The strong 
coupling constant is fixed to as{Mz) = 0.118 at the Z 
boson mass and uses the 2-loop formula for the evo- 
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FIG. 37: (color online) MC closure test of the method used to 
extract the inclusive jet pr cross section for the jet \y\ < 0.4 
bin. The full analysis was repeated treating MC events as 
data and comparing the result to the input cross section. 
Good agreement is found within the statistical uncertainties 
of fits to jet energy scale and resolution, and unfolding present 
in MC (shaded band), which are much smaller than the sys- 
tematic uncertainties in data. 



lution of as- The Pythia cross section is reweighted in 
s so that the Pythia parton shower prediction agrees 
with NLO pQCD. The correction factors for hadroniza- 
tion and the underlying event are shown in Fig. 1381 As 
shown in Fig. I36[ the measurement is well described by 
NLO QCD over eight orders of magnitude in the six y 
bins. 

To check more precisely how well the measurement is 
described by the NLO QCD theory, we display the ra- 
tio of data over theory in Fig. [39l where the theory is 
calculated using the CTEQ6.5M PDF parameterization. 
The PDF uncertainties represented as dashed lines are 
calculated using the set of 20 eigenvectors provided by 
the CTEQ Collaboration for the CTEQ6.5M PDF fits. 
Data and theory agree within experimental and theoret- 
ical uncertainties, but data seems to favor the lower end 
of the CTEQ6.5M PDF 90% confidence level uncertainty 
band. Data are also compared to the NLO QCD calcu- 
lations using the MRST2004 PDF parameterization and 
our agreement in shape is good. The experimental uncer- 
tainties are smaller than the present PDF uncertainties, 
so these data further constrain the PDFs. 

Some recent parameterizations have already used our 
measured jet cross sections described here to further con- 
strain the PDFs. As an example, we display in Fig.l^the 
ratio data over theory, where the NLO theory is calcu- 
lated using the MSTW08 NLO PDF [3l| which displays 
good agreement between our measurement and this pa- 
rameterization, with a tendency to be slightly different 
at high jet pt where the uncertainties are larger. For 
reference, we also display in Fig. UT] the ratio of data 
over theory where the theory uses the recent HERA- 
PDFvl.O PDF [32], which uses only HERA data to con- 
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strain PDFs. We notice some discrepancies between our 
measurement and the HERAPDFvl.O PDF at medium 
jet pt especially in the central region, and at high pT in 
the forward region. We also compare our data with the 
ABKM09NLO [sj parameterization in Fig. |42] and we 
notice some disagreement between our data and the pre- 
dictions in particular on the normalization. This shows 
the capability of our data to constrain further the PDFs. 
Furthermore, we compare our measurements to the re- 
cent CTIO [3J] parameterization in Fig. |43l There is 
a good agreement with data with the tendency of the 
CTIO parameterization to be higher at large pt in all 
\y\ bins. Finally, we compare our measurement with the 
predictions from the NNPDFv2.1 [sl] parameterization 
in Fig. |44] and again good agreement is found with our 
data. 

The details of the uncertainties on the inclusive jet pT 
cross section are given in Fig. 1451 The dominant un- 
certainty is due to the systematic uncertainties on the 
jet energy scale, but the unfolding and the uncertainties 
related to the resolution in jet px are also important, 
especially at high pT and high \y\. The 6.1% luminos- 
ity uncertainty is the second largest uncertainty at low 
Pt and the third largest at high pr, and leads to signifi- 
cant uncertainty in the overall normalization of the cross 
section. For a jet pr ~ 150 GeV, it is similar to the 
jet energy scale uncertainty. The uncertainties related to 
efficiencies are small everywhere. 



XIII. UNCERTAINTY CORRELATIONS 

Correlations between systematic uncertainties are 
studied in detail to increase the value of these data in fu- 
ture fits to model parton distributions and their impact 
on LHC physics predictions in particular. In total, there 
are 91 independent sources of systematic uncertainty, and 
in this section we describe the method we use to group 
those with similar impact on the shape of the cross sec- 
tion to find the principal components of the uncertainty 
without significantly impacting the overall quality of the 
data. Many of the systematic sources we describe above 
are small in magnitude and highly correlated in shape 
with other sources. 

The traditional interpretation of uncertainties to be 
independent requires that at each point the sum of all 
sources in quadrature must equal the total systematic 
uncertainty. In practice, adding in quadrature sources 
with similar shapes whose orthogonal components (de- 
fined later) are small will lose very little information com- 
pared to the full information given in the 91 different 
systematic uncertainties. 

We combine uncertainties that are correlated and of 
similar shape to reduce the number of components in the 
covariance matrix. We develop a robust systematic ap- 
proach for regrouping the sources based on the notions 
of source size, shape similarity, and orthogonality. The 
natural measure for the size of a source is the impact it 



has on the overall in the fit with the ansatz function 
when shifted by one standard deviation around the mini- 
mum. To assess the similarity in shape between different 
systematic uncertainties, we define the inner product for 
sources h and g as 



{h-9) = 



hi ■ Oi 



iGbins 



(32) 



where hi and gi are the values of two systematic uncer- 
tainties and the sum is over the pt and \y\ bins. The 
size, or magnitude, of a source h can be written using 
this notation as 



(33) 



The shape similarity of two sources h and g can be quan- 
tified by calculating their correlation, which is written in 
the notation of Eq. [32] as 



{h-g) 
\h\\-\\g 



(34) 



which varies between -1.0 and 1.0. When p =1.0, the 
sources are fully correlated, —1.0 fully anti-correlated and 
0.0 completely uncorrelated. The source g can be broken 
into a component that is fully correlated with source h 
and another component that is fully uncorrelated by con- 
sidering a linear transformation 



i — g ^ ah. 



(35) 



When the orthogonality of h and g' is defined in terms 
of the inner product. 



h _L 



iji ■ g') = 0, 



Eq. ESI and Eq. [ 



] together yield 
_ {h-9) 



(36) 



(37) 



defining g' as the orthogonal component that is fully un- 
correlated with source h. The value g' has the property 
{g' ■ g') < (g-g), {g' ■ g') ^ (g-g) is equivalent to h be- 
ing orthogonal to g, and {g' ■ g') = to h being parallel 
to g. Small values of \\g'\\ indicate that the sources can 
be combined with little impact on the freedom of the fit 
to the ansatz. 

The sources due to statistical uncertainties in fits are 
first assigned as uncorrelated. The remaining sources are 
sorted by size and are then iteratively recombined with 
other sources most similar in shape and having the small- 
est orthogonal components. The sources are combined 
when their correlation is greater than about 85% and the 
orthogonal components have a magnitude smaller than 
10% of the largest individual systematic Cmax- At the end 
of the iterative procedure, the remaining set of sources 
no longer has any pairings with an orthogonal component 
less than O.lemax- The smallest remaining sources with 
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FIG. 38: (color online) Hadronization (dashed line) and underlying event (dotted line) corrections for inclusive jet cross section 
and the product of both corrections (solid line). The uncertainty on the theory is estimated as 50% of the individual corrections 
added in quadrature. 
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FIG. 39: (color online) Measured data divided by theory for the inclusive jet cross section as a function of jet pr in the six 
\y\ bins. The data systematic uncertainties are displayed by the shaded band. NLO pQCD calculations, with renormalization 
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the data. The CTEQ6.5 PDF uncertainties are shown as dashed lines and the predictions with MRST2004 PDFs as dotted 
lines. The theoretical uncertainty, determined by changing the renormalization and factorization scales between pt/2 and 2pT, 
is shown at the bottom of each figure. 
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FIG. 40: (color online) Measured data divided by theory for the inclusive jet cross section as a function of jet pr in the six 
|y| bins using the MSTW2008 parameterization. The data systematic uncertainties are displayed by the shaded band. The 
theoretical uncertainty, determined by changing the renormalization and factorization scales between pt/2 and 2pT, is shown 
at the bottom of each figure. 
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FIG. 41: (color online) Measured data divided by theory for the inclusive jet cross section as a function of jet pt in the six 
\y\ bins using the HERAPDFl.O parameterization. The data systematic uncertainties are displayed by the shaded band. The 
theoretical uncertainty, determined by changing the renormalization and factorization scales between pt/2 and 2pT, is shown 
at the bottom of each figure. 
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FIG. 42: (color online) Measured data divided by theory for the inclusive jet cross section as a function of jet pr in the six 
\y\ bins using the ABKM09 parameterization. The data systematic uncertainties are displayed by the shaded band. The 
theoretical uncertainty, determined by changing the renormalization and factorization scales between pt/2 and 2pT, is shown 
at the bottom of each figure. 
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FIG. 43: (color online) Measured data divided by theory for the inclusive jet cross section as a function of jet pt in the six \y\ 
bins using the CTIO parameterization. The data systematic uncertainties are displayed by the shaded band. The theoretical 
uncertainty, determined by changing the renormalization and factorization scales between pt/2 and 2pT, is shown at the bottom 
of each figure. 



magnitude less than O.lemax are added in quadrature to 
the uncorrelated uncertainty. The final reduced set of 
uncertainties has 23 correlated sources (principal compo- 
nents) and one fully uncorrelated uncertainty, which is a 
significant reduction compared to the original 91 sources. 
The reduced set of 23 correlated sources and the total un- 
correlated uncertainty are provided in Ref. ^oj . 



The five leading sources from the reduced set of com- 
bined systematic uncertainties, the total uncorrelated un- 
certainty, and the total uncertainty are shown in Figs. l46l 
and ITtI in the six \y\ bins. These sources summarize the 
leading systematic uncertainties for the measurement. 
The EM scale uncertainty comes from the calibration of 
the EM calorimeter using Z — >■ e+e~ events. The photon 
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FIG. 44: (color online) Measured data divided by theory for the inclusive jet cross section as a function of jet pr in the six 
\y\ bins using the NNPDFv2.1 parameterization. The data systematic uncertainties are displayed by the shaded band. The 
theoretical uncertainty, determined by changing the renormalization and factorization scales between pt/2 and 2pT, is shown 
at the bottom of each figure. 



energy scale includes the uncertainty in the MC descrip- 
tion of the difference in the electron and photon responses 
and the uncertainty in the amount of passive material in 
front of the calorimeter, which affects the response differ- 
ence as a function of photon p^- The uncertainty in the 
high pt extrapolation is due to differences in fragmenta- 
tion models of Pythia and Herwig, which lead to an 
additional uncertainty in the high px extrapolation of 
the central response. The rapidity-intercalibration un- 
certainty summarizes the uncertainty in the relative re- 
sponse calibration between calorimeter regions. The de- 
tector showering uncertainty includes the uncertainties 
on showering, but also additional significant contribu- 
tions from other uncertainties such as sample purity and 
the difference between alternate tunes of Pythia (tunes 
A and QW). 

XIV. CONCLUSION 

In this paper, we described the measurement of the 
inclusive jet cross section by the DO experiment. The 
measured inclusive jet cross section corrected for experi- 
mental effects to the particle level in pp collisions at ^/s — 
1.96 TeV with C = 0.70 fb~^ was presented for six \y \ bins 
as a function of jet pr- The precision reached in this 



measurement is unprecedented for results from a hadron 
collider, particularly for processes dependent on gluons 
at high-x. The measurement was found to be in good 
agreement with NLO QCD calculations with CTEQ6.5M 
and MRST2004 PDFs. These resuhs will also be useful 
for any experiment at a hadron collider such as the LHC 
where the same techniques can be used to extract the 
jet energy scale with high precision and to measure the 
inclusive jet cross section. In addition, a full analysis 
of correlations between sources of systematic uncertainty 
was performed, demonstrating a useful method to reduce 
the complexities of describing numerous sources of uncer- 
tainties in the cross section, and increasing the potential 
impact of these data in global PDF fits. 
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FIG. 45: (color online) Different components of the systematic uncertainty as a function of jet pr in the six y 
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FIG. 46: (color online) Correlated uncertainties for all central regions and the ICR as a function of jet pr for four \y\ bins, 
\y\ < 0.4, 0.4 < jj/j < 0.8, 0.8 < \y\ < 1.2, and 1.2 < \y\ < 1.6 . The five largest systematic uncertainties are shown together 
with uncorrelated and total uncertainties, computed as the sum in quadrature of all sources. 
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FIG. 47; (color online) Correlated uncertainties for 1.6 < 
\y\ < 2.0 and 2.0 < \y\ < 2.4 as a function of jet pr- The 
five largest systematic uncertainties are shown together with 
uncorrelated and total uncertainties, computed as the sum in 
quadrature of all sources. 
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